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TEMPORAL EVOLUTION OF THERMAL EMISSION FROM RELATIVISTICALLY EXPANDING PLASMA 
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ABSTRACT 

Propagation of photons in relativistically expanding plasma outflows, ejected from a progenitor 
characterized by steady Lorentz factor P is considered. Photons that are injected in regions of high 
optical depth are advected with the flow until they escape at the photosphere. Below the photosphere, 
the photons are coupled to the plasma via Compton scattering with the electrons. I show here, that 
as a result of the slight misalignment of the scattering electrons velocity vectors, the (local) comoving 
photon energy decreases with radius as e' (r) oc r~'^/^ . This mechanism dominates the photon cooling 
in scenarios of faster adiabatic cooling of the electrons. I then show that the photospheric radius of 
a relativistically expanding plasma wind strongly depends on the angle to the line of sight, 6. For 
9 < r~^, the photospheric radius rph is 0- independent, while for 9 > P^^, rpfi{9) oc 9^. I show 
that the ^-dependence of the photosphere implies that for flow parameters characterizing gamma-ray 
bursts (GRBs), thermal photons originating from below the photosphere can be observed up to tens of 
seconds following the inner engine activity decay. I calculate the probability density function P{r, 9) 
of a thermal photon to escape the plasma at radius r and angle 9. Using this function, I show that 
following the termination of the internal photon injection mechanism, the thermal flux decreases as 
F^g{t) oc and that the decay of the photon energy with radius results in a power law decay of 
the observed temperature, T°^ {t) oc , with j3 = 2/3 at early times, which changes to /3 ~ 1/2 
later. Detailed numerical results are in very good agreement with the analytical predictions. I discuss 
the consequences of this temporal behavior in view of the recent evidence for a thermal emission 
component observed during the prompt emission phase of gamma-ray bursts. 

Subject headings: gamma rays:theory — plasmas — radiation mechanisms: thermal — radiative transfer — 
scattering — X-rays:bursts 



1. INTRODUCTION 

Evidence for relativistic expansion in plasma winds ex- 
ist in various astronomical o bjects, such as microquasar s 
()Mirabel fc Rodriguez Ill994l : llliellmin g fc Rupcn 1995V 
activ e galactic nuclei (AGNs; Lind & Blandford 
[19851 iGopal-Krishna et 0/71120061) and gamm a-ray bursts 
(GRBs: iPaczvnski Ill986t iGoodman |[l98l . In many of 
these objects, the density at the base of the flow is suffi- 
ciently high, so that the optical depth to Thomson scat- 
tering by the baryon-relatcd electrons exceeds unity. If 
the optical depth to scattering is high enough, the emerg- 
ing spectrum of photons emitted by radiative processes 
occurring at or near the base of the flow is inevitably 
thermal or quasi-thermal (a Wien spectrum could also 
emerge if the number of photons is conserved by the ra- 
diative processes). These photons escape the flow once 
they decouple fro m the plasma, at the photosphere (e.g., 
IPaczvnski |]l990D . 

The photosphere is usually defined as a surface in 
space which fulfills the following requirement: the optical 
depth to scattering a photon originating from a point on 
this surface and reaching the observer is equal to unity. 
Therefore, calculation of the position of this surface re- 
quires knowledge of the density profile between this sur- 
face and the observer (and, in principle, knowledge of 
the photon energy and the velocity profile, since the cross 
section is energy dependent; however, I will neglect these 
effects). Calculation of the photospheric surface in the 
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case of steady, spherically sym metric, relativistic wind 
was carried out by Abramowicz et (iL \ (|1991). In this 
work, it was found that the position of the photosphere 
has a complicated, non-trivial shape which strongly de- 
pends on the viewing angle and the wind Lorentz fac- 
tor P. This shape can be described analytically (see 
I Abramowicz et aFlliggiL eq. 3.4). For spherically sym- 
metric wind, the photospheric surface is symmetric with 
respect to rotation around the axis to the line of sight. 
Thus, I will use the term "photospheric radius" from here 
on to describe its position in space, noting that the pho- 
tospheric radius is a function of the angle to the line of 
sight, Tph = rph{9). 

While the optical depth to scattering from the pho- 
tospheric radius rph{9) to the observer is by definition 
'''{fph[d\) = 1, in fact photons have a finite probability of 
being scattered at any point in space in which electrons 
exist. Since in every scattering event a photon changes 
its propagation direction and its energy, the observed flux 
and temperature of the thermal photons depend on the 
last scattering position, scattering time, the comoving 
temperature at this position, and last scattering angle. 
A full description of the last scattering position and scat- 
tering angle can only be done in terms of probability den- 
sity function P{r,9). The probability density function is 
an extension of the standard use of the photospheric ra- 
dius as a surface in space from which thermal photons 
emerge, to consider the finite probability of a photon to 
emerge from an arbitrary radius r and arbitrary angle 9. 

Once the probability of a thermal photon to emerge at 
time t from radius r and angle 9 is known, the observed 
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flux of tlie thermal photons can be calculated. The first 
observed (thermal) photon originates from the radial axis 
towards the observer (on the line of sight). At later times 
an observer sees photons that originate from increasingly 
higher angles to the line of sight and from larger radii. 
The observed thermal flux thus varies with time. 

The observed temperature of thermal photons emerg- 
ing from radius r at angle to the line of sight is blue 
shifted due to the Doppler effect, T°^- = VT'{r). Here, 
T'(r) is the photon temperature in the comoving frame, 
T°''- is the observed temperature, V = [r(l — /?/i)]~^ 
is the Doppler factor, (3c is the fluid velocity and ^ = 
cos(6'). Photons emitted on the line of sight are blue 
shifted by Po = 2^(0 = 0) ~ 2r (the last equality holds 
for r ^ 1), while photons that originate from 6 > are 
blue shifted by Doppler factor V{9 > 0) < Vq. Note that 
I'(0) is a monotonically decreasing function of 6. 

The photon comoving temperature T'(r) varies with 
the radius below the photosphere. At r < Vph, photons 
are coupled to the flow by multiple Compton scatter- 
ing. Therefore, once thermalize, the photons comoving 
temperature is equal to the electrons comoving temper- 
ature T^i{r) (Tgj is measured in units of nieC^)- The 
electrons comoving temperature changes as they prop- 
agate downstream, due to adiabatic energy losses and 
possible internal heating mechanisms that are coupled 
to the flow. The photon temperature is thus expected to 
trace the electrons temperature. However, as I will show 
below, an additional mechanism determines the photon 
temperature below the photosphere. This mechanism is 
based on photon energy losses due to the misalignment 
of the scattering electrons velocity vectors in regions of 
high optical depth, and leads to photon comoving tem- 
perature decay as a power law in radius, T'(r) cx r^^/^ 
(in the limit of relativistic outflows, characterized by 
Lorentz factor F 1). Therefore, as long as the elec- 
trons comoving temperature does not drop faster than 
Tg;(r) oc r~^/^, the photon comoving temperature traces 
the electrons comoving temperature. However, if adi- 
abatic energy losses cause the electrons comoving tem- 
perature to decay faster than r~^/^, than the mechanism 
described below limits the photon comoving temperature 
to decay as r~^/'^. In this case, the photon comoving 
temperature does not follow the electrons temperature. 

The mechanism by which photons lose their energy 
below the photosphere is solely based on 2- and 3- 
dimensional scattering geometry. In every scattering 
event, the direction of the photon propagation vector 
slightly changes, as the average scattering angle (9) ^ 
r^^. Therefore, in every consequent scattering, a pho- 
ton is being scattered by electrons whos direction vectors 
are slightly misaligned. This misalignment, in turn, in- 
herently leads to photon energy losses, regardless of the 
comoving temperature of the electrons. This effect can 
best be understood if one considers the rest frame of the 
first scatterer. Assuming that in this frame the comov- 
ing photon energy is very low, e' ^ ra^c^ , the (comoving) 
outgoing photon energy is nearly equal to its energy be- 
fore the scattering, independent on the scattering angle. 
The misalignment of the electrons velocity vectors im- 
plies that for scattering angles which are different than 
(0, tt), the velocity vector of the next scatterer points out- 
ward. Thus, the comoving energy of the photon in the 



rest frame of the second scatterer is slightly lower than 
its comoving energy in the rest frame of the first scat- 
terer. As I will show below, for constant outfiow velocity 
characterized by Lorentz factor F ^ 1, this effect results 
in a power law decay of the comoving photon energy, 
e'(r) cx r~^/^. This decay, in turn, results in a decrease 
of the observed temperature at late times. 

Transient physical sources have a finite emission dura- 
tion, during which their inner engine is active. Therefore, 
at any given observed time t°^^ > 0, an observer sees si- 
multaneously photons originating from a range of radii 
and a range of angles to the line of sight , ^min < ^' < ^'max ■ 
The dependence of the photon comoving temperature 
on the photospheric (last scattering event) radius and 
the dependence of both the photospheric radius and the 
Doppler factor on the angle 9 to the line of sight thus 
lead to the conclusion that thermal emission (Planck 
spectrum) in the comoving frame is observed as a mod- 
ified black body. However, the observed spectrum is a 
convolution of black body spectra observed with differ- 
ent fluxes and temperatures, and as such does not differ 
much from a black body. 

As long as the inner engine is active, the observed flux 
and temperature are dominated by photons emitted on 
axis (on the line of sight). Following the decay of the 
inner engine, the flux becomes dominated by photons 
emitted from increasingly higher angles and radii (the 
curvature effect, also k nown as "the high latit ude emis- 
sion" effect; see, e.g., iFenimore. Ma.dras fc N avakshini 
1199^; Kumar fc P a.naitescu I l200d : iRvdeTPetrosia^ 
i2_002; Dermcr 20041, for discussions on this effect in op- 
tically thin emission models). Considering the optically 
thin cases, the results obtained in these works are not 
valid for the scenario considered here, of thermal emis- 
sion from optically thick expanding plasma. 

The main goal of this paper is to provide a theoreti- 
cal framework for the analysis of thermal emission from 
astrophysical transient characterized by relativistic out- 
flows, and to calculate the expected observed thermal 
flux and temperature at late times, following the decay 
of the inner engine. A key motivation to this work are 
the results obtained bv lRvde I (|2004ll2005h . of a decaying 
thermal component observed during the prompt emission 
phase of GRBs. In these works, it was shown that after 
~ 1 — 3 s the temperature of the thermal component de- 
creases as a power law in time, T"**' cx with power 
law index a ~ 0.6 — 1.1. An additional analysis (F. Ryde 
& A. Pe'er 2008, in preparation) shows that after a short 
rise, the flux of the black body component of these bursts 
also decreases with time as F'^^ cx with power law 
index (3 « 2.0 — 2.5. I show here that these results are 
naturally obtained in a model which considers the full 
spatial scattering positions and photon scattering angles 
(given by the probability distribution function P[r,9]) 
and takes into consideration the comoving energy losses 
of the photons below the photosphere, due to the slight 
misalignment of the scattering electrons. 

This paper is organized as follows. I first calculate in 
[J2] the optical depth at every point in space for scat- 
tering into angle 0, under the assumption of relativistic, 
steady outflow. From this calculation, I flnd the angular 
dependence of the photospheric radius, rph{9). The cal- 
culation in this section closely follows the treatment by 
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lAbramowicz et aUl ()1991h . although it is somewhat more 
general since in this work the results were not presented 
as a function of the viewing angle, 6. For parameters 
characterizing GRBs, it is shown in i )2.1l that thermal 
emission can be observed up to tens of seconds. The 
results obtained in this section are used in f|3] in which 
the theory of photon energy loss due to misalignment of 
the scattering direction vectors below the photosphere is 
developed. The probability density function P{r, 9) of a 
photon to be scattered from radius r and angle 9 is intro- 
duced in fj4l In this section, I calculate the expected flux 
and temperature of the thermal emission, and show that 
following the termination of the inner engine, the ther- 
mal flux and temperature decay as power law in time, 
F cx f-", and T"''- cx t"", with a = 2 and /3 = 2/3 at 
early times, modified into (3 ~ 0.5 later. The numerical 
model is presented in 33 The Monte-Carlo simulation 
provides full calculation of photon propagation in rel- 
ativistically expanding plasma, and is used to validate 
the approximations of the analytical calculations. The 
numerical results are compared to the analytical predic- 
tions, and are found to be in very good agreement. I 
summarize the results and discuss their consequences in 

m 

2. OPTICAL DEPTH AND PHOTOSPHERIC RADIUS IN 
RELATIVISTICALLY EXPANDING PLASMA WIND 



Following the treatment of lAbramowicz et al. I (|1991l ). 
I consider the ejection of a spherically symmetric plasma 
wind from a progenitor characterized by constant mass 
loss rate M, that expands with time independent veloc- 
ity V = (3c. The ejection begins at i = from radius 
r — Q, thus at time t the plasma outer edge is at ra- 
dius rout{t) = Pet from the center. For constant M and 
F, at r < rout the comoving plasma density is given by 
n'{r) = M/(47rmp?;Fr2), where F = (1 - (3^)-^/^. I as- 
sume that emission of photons occurs deep inside the flow 
where the optical depth r ^ 1 as a result of unspecified 
radiative processes. The emitted photons are coupled to 
the flow (e.g., via Compton scattering), and are assumed 
to thermalize before escaping the plasma once the opti- 
cal depth becomes low enough. In the following, I will 
assume that the plasma wind occupies the entire space, 
i.e., rout{t) oo^. 

Calculation of the optical depth is done in the following 
way. Consider a fluid element that moves in a particular 
direction (in the observer frame) with constant velocity 
V. I assume that the cross section for photon scattering 
is energy independent, and is equal to Thomson cross 
section, gt- This assumption holds as long as in the 
(local) comoving frame of the fluid, the photon energy is 
low, e' ^ TTieC^. As a consequence of this assumption, 
the mean free path of photons in the plasma comoving 
frame, /' = (nVx)"^, is independent of the fluid velocity 

^ The expansion of the plasma during the photon propagation 
impHes that for photon emitted at re < rout{t), the plasma expands 
to radius ~ r^(rout[t] — re), where t is the photon emission time, 
before the photon crosses it. For parameters characterizing GRBs, 
the optical depth obtained by the full calculation converges to the 
result obtained using the approximation Tout ^ oo on an observed 
timescale of milliseconds. Full calculation implies that during this 
early expansion stage photons emitted from angle 9 > T are 
obscured; However, this calculation is omitted here, being relevant 
only for the very early stages of the expansion, and will appear 
elsewhere. 



V, as measured in the observer frame. 

Consider the propagation of photons through the 
medium in a direction which makes an angle 9 with re- 
spect to V, in the observer frame. The mean free path 
of photons as measured in this frame is I = I' — f3fi), 
where = cos(6'). Fix two points along the light path, 
with distance ds. If the fluid is at rest, /3 = 0, the optical 
depth at distance ds is dro = n'uTds. When the fluid 
velocity is nonzero, the optical depth is 
ds V 

dr = — = -jd.To = F(l — /3fi)n' ards. (1) 

With equation [1] in hand, one can calculate the op- 
tical depth for propagation of photons in spherically 
symmetric expandi n g win d. Following the treatment of 
lAbramowicz et ai] (|l99ll ). I define a cylindrical coordi- 
nate system centered at the plasma expansion center, and 
assume that the observer is located at plus infinity on the 
z axis. Consider a photon that propagates towards the 
observer (in the +z direction) at distance r^in from the 
z axis. Its distance from the center will be denoted as r, 
r — (r^-ji„ -I- z^y^^. The photon is assumed to be emit- 
ted at point z,nin along the z axis. The optical depth 
measured along a ray traveling in the +z direction and 
reaching the observer is given by 

r(rinin,2;min) = J^^,^n' {r)aTT[l - P coa{9)]dz 

(2) 

-1/21 



■| — tan 



-(a 



1 



where cos{9) = z/r 
was used. Here, 



z/(, 



Rd = 



min + •2^)^^^! and equation [T] 
Mctt 



Equation |2l is identical to the result obtained by 
lAbramowicz et al. I ()199lD . 

At the photon emission location, the angle 0(rmin, z = 
Zmin) is the angle to the line of sight. This allows to write 
equation [2] in a simpler form. Using tan(0) = r-^i-a/ z-^i-a^, 
one obtains 

Rd ( 1 



Trr 



sin(6l) 



-P 



2^ F^ + ^l- 



The last equality holds for F 1 and small angle to 
the line of sight, 9 <C 7r/2, which allows the expansion 
sin(0) ~ 61 - 6*3/6. 

The photospheric radius is obtained by setting 
T{rph,9) = 1, 

I thus find that at small viewing angle 9 ^ the pho- 

tospheric radius is angle independent, Vph — i?d/27rF^, 
while for large angles 9 ^ F~^, the photospheric ra- 
dius is rph{9) ~ Rd9^ / Qt: .'^ The full calculation of the 
photospheric radius from equation IH as well as the ap- 
proximated solution in equation [5] are presented in figure 
[TJ It is clear from the figure that the approximate solu- 
tion is nearly identical to the exact solution in the range 
9 <1 radian. 

"* Note that often in the literature, one dimensional calculation 
is used for the photospheric radius. The result obtained, 
is correct in the limit 9 ■C r~^. 
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10 10" 
Angle to the line of sight, 9 

Fig. 1. — Normalized photospheric radius Tpf^id) / as a func- 
tion of the angle to the Une of sight, 9. The thick lines show 
the exact solution of equation |4] for F = 100 (lower line) and 
F = 5 (upper line). The dotted lines show the approximate so- 
lution Tpii/R^ = {0^/3 + r~^)/2n. The approximate result in 
equation [5] is nearly identical to the exact solution for 8 < 1. 



2.1. Characteristic time scale for observation of thermal 
emission in GRBs 

The strong angular dependence of the photospheric ra- 
dius (equation imphes that thermal photons originat- 
ing from high angles to the line of sight can be observed 
on a very long time scale following the decay of the in- 
ner engine that produces the thermal emission. Below 
the photosphere, the photons are coupled to the flow, 
therefore their velocity component in the direction of 
the flow is ~ /3c. ^ Assuming that a photon is emitted 
at < = 0, r = 0, it emerges from the photosphere at 
time t = rph{0)/(3c. Photons that propagate towards 
the observer at angle to the line of sight 9 are thus ob- 
served at a time delay compared to a hypothetical pho- 
ton that was emitted at t = 0, r = and did not suf- 
fer any time delay ("trigger" photon), which is given by 
^t°^ {e) = [rph{.e)/l3c] X [l-/3cos(0)]. 

For relativistic outflows, F 3> 1, photons emitted on 
the line of sight [9 = 0) are thus seen at a time delay 
with respect to the trigger photon. 



Ai°''(6l = 0) 



Rd 



47rF4/3c 



52i 2 S. 



(6) 



Here, M = L/Tc^ , and typical parameters characterizing 



GRBs, L = 10^^X52 ergs-i and F = IOOF2 were used. 

Photons emitted from high angles to the line of sight, 
9 3> F~^ (and 6* <C 1) are observed at a much longer time 
delay, 



t' 



ob. 



>F" 



Rd f d 



37r/3c 



(7) 



where 9 = 0.16'_i. The time scale derived on the right 
hand side of equation [7] is based on the estimate of the 
jet opening a ngle in GRB outflow, 9 < 9j ~ 0.1 (e.g., 
iBerger et alT 2003). One can thus conclude, that in rel- 
ativistically expanding wind with parameters character- 

^ One way to obtain this result is by noting that the average 
photon scattering angle (in the observer frame) is (9) ~ F"-*^. Thus, 
the photon velocity component in the flow direction is c X cos 9 ~ 
c(l — 6^/2) ~ f3c. This effect will be further discussed in i[3l below. 



izing emission from GRBs, thermal emission can be ob- 
served up to tens of seconds following the decay of the 
inner engine. 

3. PHOTON ENERGY LOSS DUE TO THE SLIGHT 
MISALIGNMENT OF THE SCATTERING ELECTRONS 
VELOCITY VECTORS BELOW THE PHOTOSPHERE 

Below the photosphere, photons undergo repeated 
Compton scattering with the electrons in the flow. For a 
jet with flnite, constant opening angle, the velocity vec- 
tors of the electrons propagating inside the jet are slightly 
misaligned. This, in turn, leads to photon energy loss via 
repeated Compton scattering as the photons propagate 
downstream. This mechanism is independent on the adi- 
abatic energy losses of the electrons, and thus dominates 
if the electrons comoving temperature decreases faster 
than r~^/'^ (see below). 

In order to calculate the photon energy loss, I as- 
sume that the comoving electrons temperature can be 
neglected, and take the limit T^^ = 0. Consider a sin- 
gle scattering event between a photon and an electron 
inside the flow. I explicitly assume that the photon en- 
ergy in the (local) comoving frame is low, e' <§; rrieC^. 
In the scattering event, the photon is being scattered 
to angle 6*', in the (local) comoving frame. Denoting 
by e' the photon energy before the scattering (the in- 
coming photon energy), its outgoing energy is given by 

= e' {l + {e'/meC^)x[l-cos{9')]]~^ ~ e'. Thus, 
the photon local comoving energy is not changed by a 
single scattering event^. 

Consider two consequent scattering events, the first 
of which occurs at radius and the second at radius 
r^+i. Due to the symmetry in the scattering direction, 
on the average the photon propagation direction is par- 
allel to the flow. Assuming that before the first scat- 
tering the photon propagation direction is parallel to 
the velocity vector of the first scatterer electron, follow- 
ing the first scattering event, the photon propagation 
direction is at angle 9^ with respect to the first elec- 
tron velocity vector (in the lab frame), where cos(6'-y) = 
[[3 + cos(6'')]/[l + l3cos{9')\. After being scattered, the 
photon travels a distance Ar^ until it is being scattered 
again. The velocity vector of the consequent scatterer 
electron makes an angle 9f,i with the velocity vector of 
the first electron (see figure [2]) . The two electrons as- 
sume to propagate at a similar Lorentz factor, F. As 
discussed above, the photon comoving energy in the rest 
frame of the first scatterer, e'^ is unchanged by the scat- 
tering process. However, Lorentz transformation to the 
rest frame of the consequent scatterer implies that the 
photon comoving energy in this frame is given by 



= e: 



[[t^ [l-/?2cOs(0eO] 

-F2/3 [cos(0,O - 1] 
-F/3sin(^,0 rn!S..)| }- 



(8) 



The photon scattering angle 9^ and the angle between 
the two electrons velocity vectors, 9^1 are related via 



tan(6'e;) 



Ari sin(0^) 



(9) 



Ti + Avi cos{9^) 

^ That is, Thomson limit is assumed for the scattering process 
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and the assumption Ar^ <^ 



A? 



Rd 



1 

f2 



Ar,: 



Rd 



1 



(13) 

Since, on the average (At^) = 1, one obtains 



Fig. 2. — Schematic view of two consequent scattering geometry 
in the lab (observer) frame. 



(see figure 

Equation [5] can be considerably simplified by not- 
ing that in the lab frame the photon scattering angle 
9.y ~ ^ 1, and that below the photosphere the aver- 
age distance traveled by the photon between consequent 
scattering is small, An/r-i -C 1 (see below). Thus, one 
can approximate tan(^?e/) « O^i ~ (Ar^/ri) sin(0^) and 
cos{9ei) w 1 — 0li/2. With these approximations, equa- 
tion [8] becomes 



-i+l 



+ 



sin^(6l^)[cos(e-,)-/3] 



l-/3cos(e.^) 



o ( AnA sin'' (8^) \ 
''^ \ ri ) l-/3cos(e^) J • 



(10) 



Of the three terms in the right hand side of equation 
[10] that contribute to the energy change between conse- 
quent scatterings, the last term is the dominant. Using 
sin(6'^) - = 0(r~i), one finds that 1 - Pcos{e^) = 
©(r-^), and cos(6'^) -(3 = 0{r-^). Therefore, the first 
two terms are of the order of (Ari/r.^)^, while the last 
term is of the order of (Ar^/r^). Thus, for (Ar^/r^) ^ 1, 
equation 1101 is approximated as 



1-/3 



Ar, 



r) 



l-/3cos(6i-.y) 



(11) 



In order to determine the energy loss of a photon it is 
thus required to estimate the ratio (Ar^/r^). Calculation 
of this ratio is done by transforming to the cylindrical 
coordinate system presented in ^21 in which the z axis 
is the photon propagation direction. The first of the 
two considered scattering occurs at location z^ia along 
this axis, and the proceeding scattering occurs at loca- 
tion 2,nax- Under these definitions, Ar.; = z,nax — ^min- 
The small scattering angle (in the lab frame) implies that 
9y w tan(0^) = riiiin/^min, and that the radii of the con- 
sequent scatterings can be approximated as ~ z^in, 
ri+i ~ Zmax (see figured]). 

Using these approximation in equation (4] one finds 
that the optical depth at the first scattering point is 
n ~ (i?d/27rzmin) X + (rmin/2min)V3], and at the 
consequent scattering point it is r^+i « (i?d/27rzinax) x 
[r^^ + (r,„in/zniax)^/3]. Writing the optical depth be- 
tween consequent scattering as Ar^ = — r^+i thus leads 
to 
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27rr2 



1 



Rd ( 1 



1 
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Equation [12] can be simplified using Ar^ 



(12) 



27 m + ^i) 



(14) 



Comparison with equation |3] shows that equation [TJ] jus- 
tifies the assumption Ari/ri <C 1 used so far, as long as 
the photon propagation occurs in region of high optical 
depths. Using further this assumption in equations [TT] 
and I14( the (local) photon comoving energy and radius 
after n scattering events can be approximated as 



Tie 



nri (q) 



(15) 



Here, e'j^, ri are the photon comoving energy and radius 
at the point in which the photon is introduced into the 
plasma, (w) and {a) are the average of the functions 



/3sin^(e-, 



5?(i^+«?)[l-/3™''(^)l 



(16) 



over the scattering events angles. Note that both oj and a 
are functions only of the photon scattering angle, 9^ (for 
constant fiow parameters, Rd and F, as assumed here). 

Equation [T5] implies that the comoving energy of pho- 
tons decreases with radius as 



e' {r) oc r <°> , 



(17) 



■^max '^mm 



The values of the functions (lu) and (a) are calculated in 
below. 

4. LATE TIME DECAY OF THE OBSERVED 
TEMPERATURE AND FLUX OF THE THERMAL 
EMISSION 

As long as the radiative processes that produce the 
thermal photons deep inside the fiow are active, the ob- 
served thermal radiation is dominated by photons emit- 
ted on the line of sight towards the observer. Once 
these radiative processes are terminated, the radiation 
becomes dominated by photons emitted off axis and from 
larger radii on a very short time scale (see eq. [6]). In order 
to calculate the observed thermal flux and temperature 
at late times, it is thus required to calculate the proba- 
bility of a photon to be emitted from radius r and angle 
to the line of sight 6. Since thermal photons are coupled 
to the flow below the photosphere, the emission radius 
of these photons is in fact the radius in which the last 
scattering event takes place. 

In the calculation below, I assume that photons are 
coupled to the flow via Compton scattering below the last 
scattering event radius, r. Below this radius, the photons 
velocity vector is, on the average, parallel to the direction 
of the flow. At the last scattering event, the photon is 
being scattered into angle with respect to the direction 
of the flow. Since below r the velocity component of the 
photon in the direction of the flow is « (3c, a photon 
emitted at time t — decouples from the plasma at time 
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tdcp — r/l3c. This photon is observed at a delay with 
respect to the "trigger" photon that was emitted ait — Q, 
r = and propagates towards the observer by t°^' = 
At°^- = r/(3c X [1 - 13 cos{e)] = ru/f3c. Denoting by 
T'(r) the comoving temperature of the photons at radius 
r, the observed temperature of the photons is T"''- = 
T'{r)V, where V = [r{l - Pfi)]-'^ is the Doppler factor, 
and fi = cos{9). The observed flux and temperature are 
thus functions of r and 9 (or u). 

The calculation of the photospheric radius presented in 
SjUgives, by definition, the radius above which the optical 
depth to scattering is equal to unity. Photons, however, 
have a finite probability of being scattered at any point 
in space in which electrons exist. Thus, the calculation 
presented in fj2] needs to be extended to include the fi- 
nite probability of photons to be emitted at any radius 
and into arbitrary angle. I introduce here calculation 
of this probability density function P(r, 6) under some 
simplified approximations, which allow full analytic cal- 
culation of the flux and temperature at late times. The 
results of an exact numerical calculation are presented in 
fj5] below. It is shown there that the analytical calcula- 
tions presented here are in very good agreement with the 
exact numerical results. 



4.1. Probability of photons to decouple from the plasma 
at radius r and be scattered into angle 9 

The increase of the photospheric radius with the an- 
gle to the line of sight, 9 (see eq. [5]) implies that the 
probability of a photon to be scattered at angle 9 is r- 
dependent. Nonetheless, as suggested by the numerical 
results in ij5] below, this dependence is hmited to a cut- 
off at a maximum angle from which photons are observed 
^max(^)7 and does not affect much the probability of pho- 
tons to be scattered to smaller angles, 9 < ^max- In the 
model below, I thus make a separation of variables to 
write P{r,9) = P{r) x P{9). The validity of this as- 
sumption (as well as the other approximations used) is 
verified by the numerical results presented in SJSI 

The probability of the last scattering event of a ther- 
mal photon to occur at radius r..r -\- Sr is calculated in 
the following way. The optical depth r(r, 9) for a photon 
scattered at radius r into angle 9 to reach the observer 
was calculated in equation U) In determining the proba- 
bility of a photon to be scattered from radius r..r -\- Sr, 
the angle into which the scattering occurs is of no im- 
portance . One can therefore write the dependence of 
the optical depth on the radius as r(r) oc r^-'^ (see eq. 
H]). This optical depth is the integral over the scatter- 
ing probability of a photon propagating from radius r to 
-|-oo, i.e., r(r) = {dT/dr)dr, from which it is readily 

found that {dT/dr)\r cx r~^. As the photon propagates 
from radius r to r + 6r, the optical depth in the plasma 
changes by St = {dT/dr)\rSr. Therefore, the probability 
of a photon to be scattered as it propagates from radius 
r to r + Sr is given by 



Psc. {r.-r -\- Sr) ~ 1 — e 



J. Sr 
St oc -ir. 



(18) 



Due to symmetry, the probability of a photon propagating 
on the hne of sight to be scattered into angle d, is equal to the 
probability of a photon propagating at angle 8 with the line of 
sight to be scattered into the line of sight. 



For the last scattering event to take place at r..r -I- Sr, 
it is required that the photon will not undergo any ad- 
ditional scattering before it reaches the observer. The 
probability that no additional scattering occurs from ra- 
dius r to the observer is given by exp(— T[r]). The prob- 
ability density function P{r) for the last scattering event 
to occur at radius r, is therefore written as 



P(r) 



(19) 



The function P{r) in equation is normalized, 
P{r)dr = 1. By comparison to equation [4l the pro- 
portionality constant is rp = rpii{9 = 0) Rd/iTrT^. 

The probability of a photon to be scattered into angle 
9 is calculated as follows. Since photons undergo multi- 
ple Compton scattering below the photosphere, I assume 
that before the last scattering event the photon propa- 
gation direction is parallel to the flow. In the following, 
I neglect the dipole approximation in the last scattering 
event angle, 9' . I thus assume that in the local comoving 
frame the scattering is isotropic, i.e., da/dfl' — Const. 
This approximation was checked numerically to be valid 
(see [JS]below). Since the spatial angle dil' = sm 9' d9' dcp' , 
the probability of a photon to be scattered to angle 9' 
(in the comoving frame) is dP/d9' oc sin^'. Integrating 
over the range < 6*' < tt gives the normalization factor 
1/2. Thus, the isotropic scattering approximation leads 
to P{9') = (sin6l')/2. 

Using Lorentz transformation to the observer frame, 
the probability of scattering into angle 9 with respect to 
the flow direction, which is the observed angle to the line 
of sight is 



P{9) = P{9')j 



9' d cos 9' d cos 



cos p dcost 



(20) 



(1-/3 cos e)2 ■ 

Using the definition u = 1 — f]cos9, equation [50] becomes 

1 



P{u) 



(21) 



Note that 1 — f3 < u < 1+/3, and the function P{u) in 



equation 1211 is normalized, ^^^^ P{u)du =^ 1. 

4.2. Temporal evolution of the observed flux and 
temperature at late times 

The diffusion model presented above implies that ther- 
mal photons emerging from the expanding plasma (i.e., 
last scattered) at radius r and into angle 9 are observed 
at time = ru/jSc. This assumption implies that all 
the photons are introduced into the plasma (by radiative 
processes occurring deep inside the flow) at the same in- 
stance, i.e., the photon injection function is assumed to 
be a (5-function in time. For finite photon injection func- 
tion, the result is a convolution of the (5-function calcula- 
tion presented here. For a quick termination of the inner 
engine, the (5-function approximation leads to a good de- 
scription of the temporal evolution of temperature and 
flux at late times, once the inner engine terminates. 

Following the decay of the inner engine, the observed 
flux at time t°^- is proportional to the probability of pho- 
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ton emission from radius r and into angle 9 by 



Fit"^-) = Fo j;-; P{r)dr P{u)du x 6 (t 



/3c 



[i;i(Zmin) - £^1(2:1, 



(22) 

At a given observed time t°^^ , the integration bound- 
aries are rmax = (3ct°^- /umiu = r2/3ci°^ /(l + /3), and 
''min = f3ct°^-/{l + (3). In evaluating the integral in the 
second line, I use z = ro/r, which lead to ^max = (1 + 
P)ro/{Pct°'' ) andzmin = (l-/3)ro/(/3ci°''-). These can be 
written with the use of normalized time, = ro(l — /3)/c 
as z„,ax = [{l+P)/{l-P)]itN/Pt°^) andz„,i„ = iw//3r^-. 
In the final formula, Ei{z) = e~*t~^dt is the expo- 
nential integral. 

At late times, t°^- ^ t^, the difference of the two expo- 
nential integrals can be written as Ei{zmin) — Ei{zma.x) ~ 
log(2;niax/'Zmin), which enables to write the temporal de- 
cay of the observed flux as 



Fit'">-)^Fo 



■log 



l + fi 
l-(3 



(23) 



I thus find that the thermal flux decays at late times as 

F{t°''-)o^r^~'^. 

In order to calculate the temporal change in the ob- 
served temperature, the power law index of the photons 
comoving energy decay with radius, resulting from the 
misalignment of the velocity vectors of the electrons be- 
low the photosphere needs to be specified. The average 
values of the functions (w) and (a) (eq. [TH]) are deter- 
mined with the use of the probability density function 
P{u) given in equation [211 By doing so, I assume that 
the conditions that led to the validity of equation!^ for 
the last scattering event (i.e., that before the scattering 
the photon propagation direction is parallel to the flow, 
and the neglection of the dipole approximation) hold for 
every scattering below and close to the photosphere. I 
further use the fact that F 1 and the average pho- 
ton scattering angle below the photosphere {9^) ^ 1 to 
approximate u = 1 - /3cos9^ ~ (^7/2 + l/2r^), and 
sin^^ ~ 9^. Using these approximations in equation 1161 
one obtains 



(a) 



2/3 



(24) 



The mean of the functions {u ^}, {u ^} is calculated 
using equation 1211 



2. 



(O = il(3 + /?2)^|r4. 



(25) 



Using equation [23 it is readily found that (lu) ~ 
(l/3)/3rV(i?d/2^), and (a) ~ (l/2)rV(-Rd/2^). Using 
these results in equation [17] leads to the conclusion that 
the comoving photon energy decays with radius as 



£ (r) ocr ~ r '"^ 



(26) 



The arguments leading to equation l261 in particular the 
requirement Ar^ /r^ ^ 1 are valid below the photosphere 
(see 521, eq. [H|) . Therefore, at radii much larger than 
ro a deviation from this law is expected. Comparison 
with the numerical results ([JSj figure[5|below) shows that 
indeed at large radii r S> rg the photon comoving energy 
becomes r independent. This result can be understood 
since above the photosphere the optical depth is smaller 
than unity, and, if a photon is being scattered at all then 
the number of scattering it undergoes is no more than one 
or two at most. In calculating the observed temperature 
at late times, I thus assume that the comoving photon 
energy decreases with radius as s'{r) oc r~'^^^ at r < 
Abrk X tq, and e'{r) cx r° at larger radii, where Abrk — 
few. 

The photons spectral distribution is thermal, resulting 
from thermalization processes occurring in regions deep 
inside the flow, which are characterized by very high op- 
tical depth. As the photons propagate outwards, their 
energy decreases, however the thermal spectrum is un- 
changed (see further discussion in |J6| below). At any 
given instance, an observer sees photons emitted from a 
range of radii and angles. Therefore, even if the comoving 
energy spectrum of the photons is thermal (black body), 
the observed spectrum deviates from black body, and is 
a grey body. Nonetheless, being a convolution of black 
body spectra, the observed spectrum is not expected to 
deviate much from black body spectra. The effective 
temperature of the observed spectra is 



ob. / j-ob. \ 



P{r)dr P{u)duT'"' {r, u)S 



ru 
Sic, 



P{r)dr P(u)du5 (t""- = ^) 

(27) 



where T°''-{r,u) = T'(r)P = T'{r)/Tu, and 



T'{r) 



r < AbrkTo, 



r>Abrkro- 



(28) 



In evaluating the numerator in equation [571 one needs 
to discriminate between two cases. At early times, t°^^ < 
Abrk^N I (i, ''max < Abrki"o, and therefore the break in the 
comoving temperature occurs at radius which is outside 
the integration boundaries. At these times, equation 1271 
becomes 

rpob.(^ob..-i-ob. ^ A i \ _ Tg f2/3 J^^.^ ^3^3 e or 

1 [t ,t SAbrktN) - J^-max ii)e-(.o/'-)dr 

rp' f'^l^r^l^ 1 ^ 

r /3ct°''- V 3/ 

[P(-i.„ax)-P(-i.„,n)] 
^ Bl(^„.„)-i5l(z„ax) • 

(29) 

Here, P{a,z) = [l/r(a)]Jp e 4" ^dt is incomplete 
Gamma function®. At early enough times, t°''- ^ tj^, 
-Zmax ^ 2niin 3> 1. Equation [29] can be put in a 
simpler form by expanding the exponential integrals 
Ei{z) c± e~^/z and the incomplete Gamma function, 
r(-l/3)P(-l/3, z) ~ -e-^z^^/^. With these approx- 

* Note that r(a) is Gamma function with the parameter a, not 
to be confused with the Lorentz factor, F. 
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imations, equation 1291 becomes 

oh. (ioh. . j-oh. ^ \ _t \ _ 



-2/3 1/3 



(30) 



At later observed times t°^^ > Abrktw/P, the break 
in the comoving temperature occurs at radius which is 
inside the integration boundaries. Sphtting the integral 
over r in the numerator of equation [57] into two, one 
obtains 



Sl(Zmi„)--El(Zmax)' 



(31) 



where 



4-2/3 
^brk 



A 



brk^ 



(32) 



Unfortunately, for the relevant time scale t°^' /t^ ^ 
lO'^, (see ^ there is no simpler analytic approximation 
to equation [311 

5. NUMERICAL CALCULATION OF THE FLUX AND 
TEMPERATURE DECAY AT LATE TIMES 

The analytical calculations presented above were 
checked with a numerical code. The code is a Monte- 
Carlo simulation, based on earlier code developed for the 
study of photon propagatio n in re lativistically expand ing 
plasma ([Pe'er fc Waxma^[2004l : [Pe'er et aL H2006U ). I 
give below a short description of the numerical code, be- 
fore presenting the numerical results and a comparison 
to the analytical approximations developed above. 

5.1. The numerical model 

I consider a three-dimensional plasma wind expand- 
ing from an initial radius that fills the entire volume 
r > ri. F ollowing the standard d ynamics of GRB out- 
flow (e.g.. [Meszaros fc Rees [[2000l ). the plasma assumed 
to accelerate up to the saturation radius rg = Fr^ , above 
which it expands at constant Lorentz factor, F. The 
plasma is assumed to be ejected at a time independent 
rate, M — L/Tc^, where L is the (observed) luminosity. 
Therefore, the plasma comoving density decreases with 
radius above rg as n'(r) = M/(4mpF/3cr^). Since be- 
low Ts, F(r) oc r, considering adiabatic energy losses the 
plasma comoving temperature at r > is given by 



FrneC^ I 4:iTr'fca J K''^^ J 

in-3 rl/4 l/6p-l/3 -2/3 



-2/3 



(33) 



Here, T^i is given in normalized units of irieC^ and the 
convention Q = lO^Qx is adopted in cgs units. 

Photons are injected into the plasma in a random 
position on the surface of a sphere at radius rj„j — 
Rd/i'^nT'^d). The depth d is taken as d = 20 in or- 
der to ensure that the probability of a photon to escape 
without being scattered is smaller than exp(— 20), i.e., 



negligible^. The initial photon propagation direction is 
random, and its (local) comoving energy at the injection 
radius is equal to the plasma comoving temperature at 
this radius. 

Given the position of the scattering event and the pho- 
ton 4- vector, the position of the next scattering event is 
calculated as follows. The code transforms the scattering 
position into the cylindrical coordinate system presented 
in Sj2l in which the scattering position is (r„iin, 2^min)- Us- 
ing 9 — tan~^(rmin/zniin), the code calculates the optical 
depth for the photon to escape, using equation[4l It then 
draws an optical depth At from a logarithmic distribu- 
tion, which represents the optical depth traveled by the 
photon until the next scattering event. If At is larger 
than the optical depth to escape, the photon assumes 
to escape, and the time difference between this photon 
and a hypothetical photon that propagated parallel to 
the last propagation direction of the photon and was not 
scattered is calculated. 

If At is smaller than the optical depth to escape, 
the next scattering position is calculated along the pho- 
ton propagation direction, i.e., it occurs in position 
(f'min, ^max)- z^ax IS Calculated such that the difference 
in optical depths between the initial scattering point and 
the next scattering point is equal to At. The scattering 
position is then transformed back to the standard Carte- 
sian coordinates. 

Given the scattering position, the electrons tempera- 
ture is drawn from a Maxwellian distribution with tem- 
perature given by equation 1331 The photon 4- vector is 
Lorentz transformed twice: first into the (local) frame of 
the bulk motion of the flow, which assumed to move at 
constant Lorentz factor F in the radial direction. Sec- 
ond, transformation into the electrons rest frame (the 
electron assumed to move randomly within the bulk mo- 
tion frame). The photon interact with the electrons via 
Compton scattering. The full Klein-Nishina cross section 
is considered in the interaction, in which the scattered 
angles of the outgoing photon are drawn. The outgoing 
photon energy and its new propagation direction are cal- 
culated in the electrons' rest frame. The program then 
Lorentz transforms the photon 4- vector back into the lab 
(Cartesian) frame, and repeats the calculation until the 
photon escapes. 

5.2. Numerical results 

The position of the last scattering events for N — 
10^ simulated photon propagation inside the expanding 
plasma is presented in figure [3] The last scattering event 
points are shown in the r — 9 plane, on top of the photo- 
spheric radius calculated in equations [4] and [5] In prepar- 
ing the plot, parameters characterizing GRBs (sec cq. 
[55]) were taken. In the figure, the last scattering event 
radius is normalized to tq = Rd/'^TrT^. Therefore, results 
obtained for arbitrary values of the free model parame- 
ters (L and F) that characterize astrophysical transients 
other than GRBs, such as AGNs or microquasars are 
similar to the ones presented. Clearly, the photospheric 
radius calculated in equation [5] gives a first order ap- 
proximation of the last scattering events radii and an- 

^ In fact, from eguations llSI and 1 2 5 1 it is found that the average 
number of scattering prior to photon escape is 2d. This result 
was confirmed numerically. 
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gles. However, it is obvious from the figure that photons 
decouple from the plasma at a range of radii and angles, 
necessitating the use of the probability density functions 
calculated in ^ 

The normalized observed flux is presented in figured! 
together with the analytical result in equation [221 In 
producing this figure, the photons are assumed to be in- 
jected as a i5-function in time. As the photons propagate 
outwards, the program traces the time in which every 
scattering event occurs. Once the photon propagation 
direction after the last scattering event is known, calcu- 
lation of the observed time is done by calculating the lag 
of the particular photon with respect to a hypothetical 
( "trigger" ) photon that propagated in a direction paral- 
lel to the final propagation direction of the photon, and 
did not undergo any scatterings. 

Physical transient sources emit during a finite time du- 
ration. Since the photon injection function is assumed 
here to be a (5 function in time, it is clear that the re- 
sults presented are valid for late time emission only, af- 
ter the central engine that produced the photon emis- 
sion had decayed. The early time rise in the flux is 
thus expected to deviate from the results presented here, 
and depend on the properties of the emission mechanism 
in physical transients. The late time power law decay, 

F{t°'' ) oc t°^' is a prediction of the model. The time 
scale in figure HI is presented in normalized units of i at , 
and thus is useful for any transient sources. 

For parameters characterizing GRBs, ~ 
10~^ Li^2^2^ ^- Since, as calculated above (see W2.1[ eq. 
[7]), thermal emission from GRBs is expected to last tens 
of seconds, the relevant time scale is t/^Ar « 10^. It is 
shown in figure [4l that on this time scale the approxima- 
tion in equations 1221 1231 arc in excellent agreement with 
the numerical results. This fact confirms the validity 
of the approximations introduced in the analytical 
calculations in 21 

The numerical results show that at very early times, 
t/tN < 10~^, the analytical result presented in equation 
[SUdoes not predict the high flux expected. This indicates 
the limitations of the model, in particular the assumption 
that the observed time depends only on two parameters, 
the last scattering event radius and angle. As the pho- 
tons diffuse below the photosphere, inevitably there is a 
small spreading in their arrival times to radius r, which 
is not considered in the analytical calculations. This dis- 
crepancy, however, is limited to very short time scales, in 
which, as discussed above, the actual nature of the inner 
engine activity determines the observed flux. 

The (local) comoving photon energy at the last scatter- 
ing event radius is presented in figure [51 It is shown that 
indeed for r < rg, e'{r) oc r~^/^, as calculated in equa- 
tion l26l This justifies the approximations that led to that 
equation. At larger radii, the approximations that led to 
equation [26l no longer hold, as the number of scattering 
a photon undergo above this radius is no more than a 
few. As a result, the photon comoving energy becomes 
r-independent. 

The decay law index of the electrons comoving tem- 
perature in equation 1331 2/3, is similar to the decay law 
index of the photon temperature calculated in fJHl In 
order to check that the mechanism described in [J3l is in- 
deed independent on the decay law index of the electrons 
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Fig. 3. — Position of the last scattering event point in 9-r plane 
for 10^ events. Tlie solid (green) line is the photospheric radius, 
calculated in equation \5\ Clearly, the last scattering events take 
place in a range of radii and angles. The photospheric radius gives 
a first order approximation to the position of these events. The 
contour lines are added to the plot in order to indicate the density 
of the emerging photons radii and angles. 

comoving temperature, additional runs with decay law 
index larger than 2/3 were performed. The results ob- 
tained in these runs were similar to the results presented 
in figure |5l 

The observed temperature as a function of time is pre- 
sented in figure |6l The numerical results are shown by 
the blue solid line, and the analytical approximation cal- 
culated in equations [29l and [31] are presented by the red 
dashed line. In preparing the plots, A^^fc = 3 was taken, 
in accordance to the numerical results of the comoving 
energy decay at large radii (see figure [H]) . Clearly, the 
analytical formula gives a good approximation to the 
numerical results, although the numerical results show 
that the power law decay of the temperature is some- 
what steeper than the analytical approximation at late 

times: T°^{t°^) oc t°^~'^ , with /3 ~ 2/3 at t/^Ar < 30 
becoming /3 « 1/2 at later times. The model presented 
here thus predicts a late time power law decay of the ob- 
served temperature, with power law index /3 that slightly 
decreases over a time scale of ~ lOOtAr from (3 ~ 2/3 to 
p ~ 1/2. For parameters characterizing GRBs, this time 
scale corresponds to few seconds. 

6. SUMMARY AND DISCUSSION 

In this paper, I addressed the question of late time 
thermal emission from optically thick, relativistically ex- 
panding plasma winds. I first showed in [21 that the pho- 
tospheric radius depends on the angle to the line of sight, 
in a non-trivial way: for 9 < F"^, rph is ^-independent, 
while for 9 > F~^, rph{9) oc 9"^ (eq. [S]). I used this re- 
sult in i i2.1l to show that for parameters characterizing 
emission from GRBs, thermal emission can be seen up 
to tens of seconds following the decay of the inner engine 
(eq. [7]). In iJ3]I showed that photons lose their energy due 
to repeated Compton scattering below the photosphere. 
The mechanism responsible for this energy loss is based 
on the geometrical effect of the misalignment between 
the velocity vectors of the electrons in the expanding 
plasma jet. This mechanism is therefore unrelated to 
other mechanisms discussed so far in the literature (e.g., 
adiabatic expansion, or the mechanism leading to Kom- 
paneets equation) . I showed that as a result of this mech- 
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Fig. 4. — Observed thermal flux as a function of the observed 
time. The solid (blue) line is the numerical simulation result, and 
the dash (red) line is the analytic approximation in equation 1221 
Time is given in units of normalized time, tjv- For t/tjv ^ 1, the 

thermal flux decays as F{t°''-) oc t°''-~^ 
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Fig. 5. — Photon (local) comoving energy as a function of radius, 
at the last scattering event radius. Numerical results are presented 
by the solid (blue) line. The results indicate that for r/rp < 1, 
s'{r) oc r~^/^, in accordance with the results of equation 1261 At 
larger radii, the arguments leading to the result in equation 1261 
no longer hold, and the photon (local) comoving energy becomes 
r-independent. 

anism, the (local) comoving temperature of the photons 
decreases below the photosphere as e'{r) oc r~^^^ (eq. 
[26|) . I introduced in []4]the probability density function 
P(r, e) = P{r) X P{e) that extend the definition of a 
photosphere to include the actual positions and angles 
in space from which thermal photons decouple from the 
plasma. Using these functions, I calculated the temporal 
decay of the observed flux (eqs. [22l [23l) and temperature 
(eqs. [29l [3T|) of the thermal emission, and showed that 
both decay as a power law in time, following the decay of 
the inner engine that produces the thermal photons. The 

flux decays as F{t°^-) on t°^- " with a and the tem- 
perature decays as T{t°^-) oc t°^-^^ with /3 = 2/3 at early 
times which later changes to /3 ~ 1/2. The analytical re- 
sults were confirmed with the results of the numerical 
simulation presented in SJSl 

The results presented here can account for the recent 
observations of thermal emission th at accompan i es lon g 
duration GRBs. As was shown by iRvde I (|2004 I2005D . 
after ~ 1 — 3 s the temperature of the thermal component 
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Fig. 6. — Observed temperature as a function of the observed 
time. The solid (blue) line is the numerical simulation results, and 
the dash (red) line is the analytic approximation in equations 1291 
and l31l with ylfj^j, = 3. The numerical results are slightly smoother 
than the analytical approximation. The temperature decays as a 
power law in time, with power law index /3 2; 2/3 at t/tjv ^ 10, 
and /3 Si 1/2 at later times. 

decreases as a power law in time, T°^- oc with power 
law index a ~ 0.6 — 1.1. An additional analysis (F. Ryde 
& A. Pe'er 2008, in preparation) shows that after a short 
rise, the flux of the black body component of these bursts 
also decreases with time as oc t^^ , with power law 
index j3 « 2.0 — 2.5. These results are thus naturally 
reproduced by the model presented here. 

A key consequence of the model presented here is that 
thermal emission at early times (before the observed 
break) is dominated by thermal photons originating from 
the photosphere on the line of sight. Therefore, obser- 
vation of the thermal emission at early times, when the 
inner engine is still active, gives a direct measurement 
of the temperature and flux of photons emitted from the 
photospheric radius on the line of sight, tq = rph{0 = 0). 
This is the innermost radius from which information can 
reach the observer. 

The interpretation presented here has a direct implica- 
tion in the study of relativistic outflows. For GRBs with 
known redshift, early time (before the break) observa- 
tion of the temperature and flux of the thermal com- 
ponent enabled direct determination of two of the least 
restricted parameters of the fireball model: the bulk mo- 
tion Lorentz factor, F and the radius at the base of the 
flow (jPe'er et al. II2007D . Being based on thermal emis- 
sion only, the method presented in this paper is insensi- 
tive to many of the inherent uncertainties in former meth- 
ods of determining the values of these parameters. Fu- 
ture measurements with the upcoming GLAST satellite 
will enable to increase the sample of GRBs with known 
redshift from which thermal emission component is iden- 
tified, to further test the model presented here and to 
gain statistics on the values of the fireball model param- 
eters. 

In addition to the prompt emission phase in GRBs, 
thermal activity may occur as part of the flaring activ- 
ity observed in the ea r ly afterglow phase of many GRBs 
(jBurrows et al. 1 120051: iFalcone et oU 120071 ) . The exact 
nature of these flares is currently not yet clear. As it 
is plausible that the flares result from renewed emission 
from the inner core, a renewed thermal emission may oc- 
cur. Analyzing this emission in a method similar to the 
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one described here and bv lPe'er eT may thus 

provide information on the flow parameters during the 
late time flaring activity. 

The relevance of the results obtained here is not lim- 
ited only to emission from GRBs, but also to emission 
from any transient phenomenon characterized by rela- 
tivistic outflow, such as AGNs and microquasars. Pro- 
vided there is a source of photons deep inside the flow, 
following the decay of this source the decay laws of the 
thermal flux and temperature derived above hold for any 
such object. An important point here, is that the nature 
of the mechanism that produces the radiation is of no 
importance, as long as it occurs deep inside the flow so 
that the photons thermalize before they escape. 

In this work, I assumed that the electrons are cold (in 
the comoving frame), and that the electrons and pho- 
tons interact only via Compton scattering. If this is 
not the case, due, e.g., to some dissipation mechanism 
that produces energetic electrons at different regions of 
the flow, than Compton scattering with energetic elec- 
trons will lead to modification of the thermal spectrum. 
This c a ,se was extensively studi ed by iRees fc Meszarosi 
er et "^ (f2005ll200l . As was shown in these 
works, the thermal photons in this case serve as seed 
photons to Compton scattering that produces high en- 
ergy, non-thermal spectrum. However, if the optical 
depth in which the energetic electrons are introduced 
into the flow is smaller than ~ unity, then the ther- 
mal component can be separated from the non-thermal 
one (jPe'er et al. I l2006l ) . Note that this is exactly the 
case in the internal collision model of GRBs: internal 
shocks can only occur at radii larger than the spreading 
radius, Tspread — ^^s — 10^^ ''i,8r| Cm, which is similar 
to ro = 6 X 10^^ L^2^2^ cm. Thus, the optical depth 
in the region where internal shocks can occur is not ex- 
pected to exceed a few. I can thus conclude that thermal 
emission is expected to be observed in GRBs under the 
assumptions of the internal collisions model. 

The late time thermal emission predicted here es- 
sentially arises from emission off the line of sight. It 
is thus similar in nature to the high latitude emission 
discussed in the literature, in the context of GRB af- 
terglow emission (iFe nimore . Madras fc Navakshin"]ll996l : 
IWoods fc Loeb1 [l999: Km nar fc Panaitescu 1120001 ). All 
these works, however, treated the optically thin case, 
which is relevant for the afterglow emission phase from 
GRBs. The work presented here differs by treating ther- 
mal emission from optically thick plasmas, characterizing 
the very early stages of emission from GRBs. 

One of the key findings in this work is the new mecha- 
nism in which photons lose their energy below the photo- 
sphere. This mechanism differs than other mechanisms 
discussed in the literature so far for radiative cooling be- 
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low the photosphere. The result obtained, e'{r) oc r^^/^ 
holds for relativistic jets characterized by constant (r- 
independent) jet opening angle. For jets in which the 
jet opening angle is r-dependent, a different power law 
decay in the photon energy is expected. 

In the calculation of the photon energy loss presented 
in 121 I neglected the electrons temperature. As the 
photons propagate downstream, their comoving temper- 
ature cannot be lower than the comoving temperature 
of the electrons. The electrons temperature decreases 
due to adiabatic expansion, which result in a decay of 
the electrons temperature as a power law in the comov- 

ing plasma volume, T^^ oc V' (for relativistic elec- 
trons) . Adopting the fireball model of GRBs (for review, 
see, e.g. lMeszarosir200Q) . above the saturation and below 
the spreading radii of the fireball, the comoving volume is 
V^' cx , resulting in a decrease of the comoving electrons 
temperature as T^i{r) oc r"^/^. Above the spreading ra- 
dius, the comoving volume increases as T^' oc r'^, which 
implies T'^i{r) oc r^^ . In any of these regimes, the elec- 
trons temperature decreases with the radius at least as 
fast as the photon temperature. 

The mechanism presented here for photon energy loss 
has some resemblance to adiabatic energy losses, as the 
photon temperature is converted into work done on the 
electrons. However, it is a different mechanism having a 
different origin. Adiabatic energy losses occur once the 
plasma expands, and its volume increases. As opposed 
to that, in the scenario considered here, the volume in 
which photons interact with the electrons (the volume 
below the photosphere) does not expand, since for con- 
stant flow parameters (M, F) the photospheric radius is 
time independent. In addition, as discussed above, the 
decay law of the photon energy is independent on the 
electrons temperature (as long as the electrons comov- 
ing temperature is not higher than the photon comoving 
temperature). The fact that between the saturation ra- 
dius and the spreading radius in GRBs the decay law 
of the photons and electrons comoving temperature is 
similar, is thus a coincidence. 
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